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ABSTRACT 

A linear-stability analysis is performed on the interface formed during the directional solidification 
of a. dilute binary alloy in the presence of a time-periodic flow. In general, the flow is generated 
by the elliptical motion of the crystal parallel to the interface. The presence of the flow can either 
stabilize or destabilize the system relative to the case without flow with the result depending on the 
frequency and amplitude of the oscillations as well as the properties of the material. In particular, 
we find that the morphological instability present in the absence of flow can be entirely suppressed 
with respect to disturbances of the same frequency as the oscillation. 

INTRODUCTION 

The manufacturing of crystals with uniform material properties is frequently hampered by the pres- 
ence of morphological instabilities during the solidification of multi-component materials. These 
nonuniformities result from an interaction between surface morphology and the concentration gra- 
dients created by solute rejection. In order to eliminate these nonuniformities, it is necessary to 
suppress the instability. In the context of directional solidification, this can be dqne by producing 
materials with sufficiently low solute concentration, or by controlling the speed of the solidification 
front 1 . 

When it is undesirable to operate within these restricted ranges of parameters, some other 
method of stabilization is necessary. In the 1960’s, Hurle suggested that flow in the melt, either 
forced or resulting from natural convection, might be used to stabilize the interface. Since then, a 
number of studies have investigated the effect of various flows on morphological stability. 

Delves 2 studied solidification into a Blasius boundary layer, using a parallel-flow approximation. 
Coriell, McFadden, Boisvert and Sekerka 3 investigated solidification into a plane Couette flow. 
Forth and Wheeler 4 and Hobbs and Metzener 5 investigated solidification into a parallel flow with 
the asymptotic suction profile. McFadden, Coreill, and Alexander 6 and Brattkus and Davis 7 looked 
at solidification into a stagnation point flow. Brattkus and Davis 8 also looked at solidification over 
a rotating disc. Merchant and Davis 9 studied solidification into a temporally modulated stagnation 
point flow, and Schulze and Davis 10 considered solidification into a compressed Stokes boundary 
layer (CSL). For a more extensive review of the role of convection in solidification see Davis 11 . 

Here, we follow up on the work studying solidification into a CSL. This flow can be generated 
by oscillating the crystal parallel to the interface. 

GOVERNING EQUATIONS 

We consider the directional solidification of a dilute binary mixture at constant speed V. We choose 
a coordinate system with x-axis located at the mean position of the crystal interface, moving with 
the front, and a z-axis that is fixed in the laboratory frame of reference. The equations governing 
the system in the fluid region are the Navier-Stokes, continuity, and solute diffusion equations: 
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To simplify the analysis, we neglect latent heat and density changes, and we assume equal densities 
and thermal properties between the two phases. We also assume that heat diffuses much faster 
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than solute, and, in this limit, the temperature field is fixed and depends linearly on the vertical 
coordinate — the frozen temperature approximation 12 . We begin by considering the problem in two 
dimensions. 

We have nondimensionalized the equations by scaling the fluid velocity with the amplitude of 
the crystal oscillations U\ the spatial coordinates with the diffusion boundary layer thickness D/V, 
where D is the solute diffusivity; time with the frequency of the crystal oscillations w; and the 
temperature and solute concentrations using the product of the diffusion boundary layer thickness 
and their respective gradients at the interface. 

The nondimensional parameters that appear in the equations and boundary conditions are the 
morphological number M, the surface energy parameter F, the Schmidt number 5, the nondimen- 
sional frequency f 2, the segregation coefficient k, and a parameter measuring the amplitude of the 
lateral oscillations e. 
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Here to is the liquidus slope in the phase diagram of the binary mixture, Coo is the far-field solute 
concentration, G is the thermal gradient, T m is the melting temperature of the solvent, 7 is the 
surface free energy, Ly is the latent heat per unit volume, and v is the kinematic viscosity. 

The interfacial conditions, evaluated at 2 = h(x, y,t) are as follows: 


u = cost (2.10) 

w = 0 (2.11) 

C = M~ 1 h-2TH (2.12) 

[1 4- Clht ecos t hx\ [1 4- (A: — 1)C] = C z ~ C x h x , (2.13) 


where H is the mean curvature of the interface. 
The basic state for this system takes the form 


u — 

e Bz cos(t — Az ) 

(2.14) 

w = 

0 

(2.15) 

c = 

1 - e~ z 

(2.16) 

h = 

0, 

(2.17) 


where A and B are known real constants that depend on the Schmidt number and forcing frequency. 

To analyse the response of this state to infinitesimal perturbations, we disturb each of these 
quantities, and separate the disturbances into normal modes of the form 


<j> = <f> 4- <j>(z, t)e iax e at 4- c.c. (2.18) 

Here we are seeking time-periodic eigenfunctions with the same period as the basic state, and o is 
the Floquet exponent. If the real part of a is not zero, then the disturbances will experience a net 
growth or decay over one period. 

Routine algebraic manipulations result in two linearized disturbance equations for the solute 
concentration and vertical component of the velocity. 
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SMALL AMPLITUDE OSCILLATIONS 

In this section we find conditions on the morphological number M such that the system is neutrally 
stable in the limit of small forcing, e — ► 0. In this limit, we find that a = 0 on the neutral curve, 
and we expand the inverse morphological number as 

M 1 — Mo -h £ 2 M 2 + . . . , (3.19) 

where Mo corresponds to the no-flow, or Mullins and Sekerka, result. 

In Figure 1 we plot the correction to the no- flow result for the CSL as a function of the 
disturbance wavenumber for various forcing frequencies. If M 2 < 0 for a given wavenumber, a 
disturbance with that wavenumber will be stabilized relative to the case with no flow. Notice that 
for Q = 1 all finite wavelengths are stabilized. The Schmidt number and segregation coefficient 
have been chosen to represent a lead- tin alloy for all of the figures in this paper. 

In Figure 2 we map out the regions of the a-Q, plane where the flow has a stabilizing influence on 
the interface (M 2 < 0). Notice that there is a range of frequencies for which all finite wavenumbers 
are stabilized. We refer to this as a window of stabilization. 

It turns out that the CSL is only able to stabilize in two dimensions. In three dimensions, only 
disturbances with wavevectors aligned with the parallel flow will be stabilized while those with 
wavevectors perpendicular to the flow are unaffected — thus the flow acts as a pattern selection 
mechanism only. 

SMALL AMPLITUDE NONPLANAR OSCILLATIONS 

Motivated by the work of Kelly and Hu 13 on Benard convection, we attempt to extend the stabilizing 
influence of the CSL to three dimensions by considering the influence of nonplanar oscillations. 
Specifically, we consider the effect of adding a second oscillation perpendicular to the first which 
has the same frequency, but may differ in amplitude and phase. This corresponds to translating 
the crystal in elliptical orbits parallel to the interface. 

For the correction due to a weak flow, we get 

M 2 — [(cos 9 + Xcos'/sin 9) 2 + X 2 sin 2 9sin 2 ^]M^ D ^ , (4.20) 

where M is the result for oscillations in one direction only (CSL), 9 is the angle that the 
disturbance wave vector makes with the x-axis, A is the amplitude ratio of the oscillations, and 7 
is the phase difference between the oscillations. 

If the two oscillations are either perfectly in or out of phase (corresponding to a phase angle 
that is an integral multiple of tt), there will be some 9 C for which the correction to the no-flow 

result, M2, is zero. This means that the flow will be able to destabilize the interface if > 0, 

but acts only as a pattern selection mechanism if M^ D ^ < 0, with disturbances oriented at the 
angle 9 C being the least stable. This is a degenerate case corresponding to a single oscillation along 
an axis lying in the x-y plane. The result is therefore equivalent to that for the CSL. 

When the phase between the oscillations, 7, is not a multiple of tt, the factor appearing in 

front of M^ D ^ is positive, and it is possible to stabilize an arbitrary three-dimensional disturbance 

provided M^ D ^ < 0. Thus the ability of the nonplanar oscillations to stabilize the interface is 
predicted by the results for the CSL, and the window of stabilization is the same for both cases. 

The maximum stabilization occurs when the phase angle is tt/ 2, with the stabilization being 
greater along whichever axis has the largest amplitude oscillations. If the amplitude ratio, A, is 
one, there will be no preferred direction for cell orientation. 

FINITE AMPLITUDE OSCILLATIONS 

In the previous section we showed that some degree of stabilization can be achieved by using small- 
amplitude nonplanar oscillations with a frequency lying within a given range. In this section we 
show that the Mullins and Sekerka interfacial mode can be entirely suppressed if these oscillations 
have sufficiently large amplitude. 
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Following the method used by Hall 14 for the Stokes layer, we find an exact solution for our 
system in the form of a Fourier series in time. We then truncate this series to produce numerical 
results. 

In Figure 3 we plot M~ l versus a for a forcing frequency fi = 10, which lies within the window 
of stabilization, and we see that by taking epsilon sufficiently large we can force the neutral curve 
below the x-axis. In Figure 4 we plot M -1 versus e for the same frequency and with the wavenumber 
fixed at the critical value for the instability in the absence of flow. Notice that when e exceeds 
about 60, the instability is completely supressed. 

CONCLUSION 

We have found that nonplanar oscillations of the crystal during directional solidification can be used 
to entirely suppress the Mullins and Sekerka instability provided the frequency of these oscillations 
is within a given range and that their amplitude is sufficiently large. We have found that this 
window of stabilization persists for a wide range of material parameters and operating conditions. 

We have searched for other modes of instability, especially subharmonics, but have not located 
any. If there are no other modes of instability, the modulation proposed here may provide a practical 
means for stabilizing the interface during directional solidification. 
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Figure 1: M 2 vs. a for S = 81.0, k = 0.3 and various ft. M 2 is independent of, T. When M 2 is 
above the x-axis, the influence of the flow is destabilizing, and when M 2 is below the x-axis, the 
influence of the flow is stabilizing. As ft — ► 00 , M 2 —*■ 0, and as ft — ► 0, M 2 approaches a steady 
state. Note that for ft = 1 the flow is stabilizing for all wavenumbers. 



Figure 2: Directional solidification into Stokes layers Regions of the a-ft plane where the 
flow stabilizes (S) or destabilizes (D) the interface relative to the case without flow. S = 81.0 and 
k = 0.3; result is independent of T. 
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Figure 3: M~ l vs. a for k = 0.3, S = 81.0, fl = 10.0, and e = {0,20,40,60}. As e is increased the 
flow stabilizes the interface, and for e « 60 the Mullins and Sekerka interfacial instability is entirely 
suppressed. 
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Figure 4: M~ l vs. e for k = 0.3, S — 81.0, 12 = 10.0, and a fixed at the critical value for the no-flow 
system. Notice that the stabilizing trend begins to reverse if the amplitude of the flow oscillations 
is too great. 
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